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We  investigate  theoretically  the  effect  of  layering  of  concrete  shields  using  a  wide  class  of  semi-empirical 
models.  We  have  found  that  (i)  the  ballistic  limit  velocity  (BLV)  of  the  multilayered  shield  does  not  depend 
on  the  order  of  the  plates  in  the  shield,  (ii)  monolithic  shield  is  superior  to  any  layered  shield  with  the 
same  thickness,  and  (iii)  the  largest  decrease  of  the  BLV  occurs  when  a  shield  is  divided  into  a  number  of 
plates  having  the  same  thickness. 
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1.  Introduction 

There  is  no  consensus  on  the  effect  of  layering  on  the  ballistic 
properties  of  shields  although  this  topic  has  attracted  considerable 
interest  among  the  researchers  in  the  field  for  a  long  time.  Most 
of  the  studies  consider  the  metallic  shields;  more  detailed  infor¬ 
mation  can  be  found  in  the  surveys  (Ben-Dor  et  al.,  2005;  Corbett 
et  al.,  1996),  monograph  (Ben-Dor  et  al.,  2006a)  and  some  recent 
publications  (Ben-Dor  et  al.,  2006b, c;  Dey  et  al.,  2007 ;  Liaghat  et  al., 
2005;  Teng  et  al.,  2007,  2008;  Zhou  and  Stronge,  2008). 

There  are  only  a  few  publications  on  this  subject  regarding  the 
concrete  shields.  It  is  noted  in  TM  5-855-1  (1986)  that  a  thick  slab 
manufactured  of  several  layers  generally  has  smaller  perforation 
resistance  than  a  single  slab  with  the  same  thickness.  Kojima  (1991 ) 
performed  experiments  with  a  spherical-nosed  projectile  having 
velocity  200  m/s  and  striking  the  reinforced  concrete  plates.  He 
tested  a  monolithic  shield  with  the  thickness  18  cm  and  double¬ 
layer  shields,  9  cm  +  9  cm  and  6  cm  + 12  cm.  The  conclusion  was  that 
the  impact  resistance  of  a  double-layer  shield  is  inferior  to  that  of 
a  monolithic  shield.  Shirai  et  al.  (1997)  investigated  experimentally 
and  numerically  the  impact  resistance  of  the  reinforced  concrete 
plates  against  the  flat-nosed  projectile  at  impact  velocity  of  about 
170  m/s.  They  compared  the  protective  effectiveness  of  the  mono¬ 
lithic  shield  with  the  thickness  of  9  cm  and  double-layer  shields 
having  different  structure,  4.5  cm +  4.5  cm  and  3  cm +  6  cm.  It  was 
found  that  double-layer  shields  had  a  higher  impact  resistance  than 
monolithic  shield.  Amde  et  al.  (1997)  conducted  computer  simu¬ 
lations  using  empirical  models  for  concrete  and  steel  plates.  They 
arrived  at  the  conclusions  that  (i)  composite  barriers  manufactured 
from  different  materials  are  more  efficient  than  layered  barriers 
manufactured  from  the  same  material  and  (ii)  multiple  barriers  are 
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less  efficient  than  single  layer  panels  having  the  same  thickness. 

Approximate  models  are  widely  used  for  estimating  protective 
properties  of  concrete  shields  (for  more  details,  see  Adeli  and  Amin, 
1985;  Bangash,  2009;  Guirgis  and  Guirguis,  2009;  Kennedy,  1976; 
Li  et  al.,  2005,  2006;  Sliter,  1980).  However,  to  the  best  of  our 
knowledge,  there  are  no  publications  where  analytical  methods  are 
applied  to  the  analysis  of  the  effect  of  layering.  The  main  goal  of  the 
present  study  is  to  close  this  gap. 

2.  Mathematical  models  of  the  layered  shield 

2.1.  Model  based  on  an  equation  for  impact  resistance 

Consider  high  speed  normal  penetration  by  a  rigid  sharp  striker 
into  a  concrete  shield.  Hereafter,  the  coordinate  h  denotes  the 
instantaneous  depth  of  penetration  (DOP)  which  is  defined  as  the 
distance  between  the  leading  edge  of  the  nose  of  the  impactor  and 
the  front  surface  of  the  shield.  According  to  an  approximate  theory 
of  penetration  which  takes  into  account  only  the  local  impact  effects 
(Kennedy,  1976),  the  instantaneous  resistance  force  is  a  function  of 
the  depth  of  penetration  h  and  the  instantaneous  impactor  velocity 
v  and  can  be  written  in  the  following  form: 

RtKv)=g{h)f{v),  (1) 

where g(/i)  and  f(v)  are  functions  determining  the  adopted  penetra¬ 
tion  model.  Particular  expressions  for  these  functions  are  suggested 
by  Kennedy  (1976),  Degen  (1980),  and  Riera  (1989). 

The  equation  of  motion  of  the  impactor,  m(d2h/dt2)  =  -R,  can  be 
rewritten  as  follows: 

mv  (fff)  = -g{h)f(v),  (2) 

where  the  instantaneous  velocity  of  the  impactor  v  is  considered  to 
be  a  function  of  h,  and  m  is  the  mass  of  the  impactor. 

The  model  determined  by  Eqs.  (1)  and  (2)  was  suggested  to 
describe  penetration  into  a  semi-infinite  shield.  In  order  to  adapt 
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Nomenclature 

b  depth  when  perforation  occurs;  thickness  of  a  shield 

A 

b  depth  of  penetration 

bW  depth  when  perforation  of  the  jth  plate  occurs; 

thickness  of  the  jth  plate 
BLV  ballistic  limit  velocity 

D  diameter  of  the  cylindrical  part  of  the  impactor 

DOP  depth  of  penetration 

e  perforation  thickness  of  a  finite  shield  normalized 

byD 

F  energy  dissipated  during  perforation 

/  function  in  the  expression  for  resistance  force  in  Eq. 

(1) 

F  function  in  Eq.  (6)  or  F(z)  =  z2 

g  function  in  the  expression  for  resistance  force  in  Eq. 

(1) 

G  function  in  Eq.  (6)  or  function  in  the  model  given  by 

Eqs.  (5)  and  (7) 

G*  function  in  Eq.  (20) 

h  instantaneous  depth  of  penetration 

H(v/X),  H elements  of  the  Hessian 
k  parameter  of  the  model  in  Eq.  (3) 

m  mass  of  the  impactor 

N  number  of  plates  in  the  layered  shield 

NDRC  National  Defense  Research  Committee 
R  resistance  force  in  Eq.  ( 1 ) 

t  time 

T  thickness  of  a  shield 

t  perforation  thickness 

Tsum  total  thickness  of  the  layers  in  the  shield 

T1^  thickness  of  the  jth  plate 

v  instantaneous  velocity  of  the  impactor 

Vjji  ballistic  limit  velocity  of  the  monolithic  shield 

vimp  impact  velocity  of  the  monolithic  shield 

vres  residual  velocity  of  the  monolithic  shield 

Vfji  ballistic  limit  velocity  of  the  layered  shield 

Vimp  impact  velocity  of  the  layered  shield 

Vres  residual  velocity  of  the  layered  shield 

V^p  impact  velocity  of  the  jth  plate 

vjll  residual  velocity  of  the  jth  plate 
x  depth  of  penetration  into  semi-infinite  shield 

z*  parameter  in  Eq.  (20) 

Greek  letters 

A  difference,  A  =  F(Vb[)  -  F(%)  in  Eq.  (29) 

A  criterion  of  optimization  in  Eq.  (33) 

A  criterion  of  optimization  in  Eq.  (35) 

li  parameter  in  Eq.  ( 10) 

&  function  determining  the  model  in  Eq.  (10) 

co  function  determining  the  model  in  Eq.  (10) 


this  model  for  a  shield  having  a  finite  thickness,  we  assume  that 
Eq.  (2)  describes  motion  of  the  impactor  until  arriving  at  the  depth 
h  =  b  where  perforation  occurs.  Note  that  because  of  scabbing,  b  is 
smaller  than  the  thickness  of  the  shield,  T.  We  assume  that 

b  =  1<T ,  (3) 

where  k<  1.  Formula  (3)  is  similar  to  the  relation  between  the  per¬ 
foration  thickness  T,  the  maximum  thickness  of  the  shield  which  is 
completely  penetrated  by  an  impactor  with  a  given  impact  velocity, 

A 

and  the  reference  depth  of  penetration  b,  the  depth  of  penetra¬ 
tion  into  a  semi-infinite  shield  by  penetrator  with  the  same  impact 


velocity.  In  particular,  according  to  the  modified  Petry’s  formula  this 
relation  is  b/T  =  0.5  (Kennedy,  1976). 

We  consider  a  range  of  impact  velocities  vimp  whereby  the  pro¬ 
jectile  perforates  the  shield.  Integrating  Eq.  (2)  by  separation  of  the 
variables  from  the  initial  conditions  at  the  impact  (h  =  0,  v  =  vimp) 
to  the  final  conditions  corresponding  to  penetrating  at  the  depth 
h  =  b  =  kT with  the  velocity  vreS}  we  obtain  the  following  relationship: 


1 Vres 


m 


7 imp 


vdv 

m 


g(h)dh. 


Eq.  (4)  can  be  rewritten  as  follows: 


F(vimp)  ~  F(Vres)  —  G(b), 


(5) 


where 


F(z)  =  m 


vdv 

WY 


G(z)  =  /  g(h)dh. 


(6) 


The  ballistic  limit  velocity  (BLV)  %  is  defined  as  the  minimum 
impact  velocity  of  the  impactor  required  for  perforating  the  shield. 
Assuming  that  perforation  occurs  when  the  impactor  emerges  from 
the  shield  with  a  zero  residual  velocity  and  taking  into  account  that 
F(0)  =  0,  Eq.  (5)  implies  that 


F(vbi)  =  G(b).  (7) 

The  shortcomings  of  the  model  determined  by  Eqs.  (1)  and  (3) 
are  associated  with  physical  interpretation  of  function  g(h)  and  the 
lack  of  validation  of  this  model  for  finite  thickness  shields. 

On  the  other  hand,  describing  the  motion  of  a  projectile  inside 
the  shield  is  redundant  in  this  study  since  what  is  needed  is  only  the 
relationship  between  the  impact  velocity  and  the  residual  velocity 
of  the  impactor.  In  the  following  we  validate  the  model  determined 
by  Eqs.  (5)  and  (7)  with  F(z)  =  z2  and  k  =  1  which  without  using  the 
resistance  model. 

2.2.  Model  based  on  the  equation  of  energy  conservation 

Based  on  the  frequently  used  assumption  (see,  e.g.,  Barr,  1990; 
Li  et  al.,  2006)  that  energy  dissipated  during  perforation,  F,  is  inde¬ 
pendent  of  the  impact  velocity,  equation  of  energy  conservation  can 
be  written  as  0.5mvfmp  -  0.5 mvjes  =  E.  Since  vimp  =  vb[  if  vres  =  0,  it 

can  be  concluded  that  F  =  OAmv^  and,  consequently,  the  follow¬ 
ing  relationship  between  the  impact,  residual  and  ballistic  limit 
velocities  is  valid: 


2  2  2 

v  —  vi0<.  =  V7,. 
imp  res  b/ 

(8) 

Eq.  (8)  coincides  with  Eq.  (5)  when 

F(z)=z2,  k  =  l(b  =  T). 

(9) 

For  studying  the  behavior  of  the  function  G(F)  =  G(fi),  we  use  well 
known  semi-empirical  models  describing  the  dependence  between 
v2h[  and  the  thickness  of  the  shield  (Adeli  and  Amin,  1985;  Li  et  al., 
2005).  After  some  algebra,  corresponding  formulas  can  be  repre¬ 
sented  in  the  following  form  that  is  convenient  for  our  analysis: 

Vj,  = /ztf'(x),  x  =  <w(e),  (10) 

where  /z  is  parameter  which  generally  depends  on  the  mechani¬ 
cal  properties  of  the  material  of  the  shield,  mass  and  diameter  of 
the  impactor,  D,  W  and  co  are  functions  determining  the  model, 
x  and  e  are  the  depth  of  penetration  into  a  semi-infinite  shield 
and  the  perforation  thickness  of  a  finite  shield  normalized  by  D, 
correspondingly. 

The  particular  expressions  for  the  functions  ^  and  co  for  a  num¬ 
ber  of  suggested  models  are  summarized  as  follows: 
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(1)  Ballistic  Research  Laboratory  (BRL)  model: 

V{x)  =  xx-5,  cu(e)  =  0.769e;  (11) 

(2)  Army  Corps  of  Engineers  (ACE)  model: 

<F(x)  =  (x  -  0.5)1'33,  a>(e)  =  0.81e-  1.06,  3  <  e  <  18; 


(12) 


(3)  Modified  National  Defense  Research  Committee  (NDRC)  model: 


nx) 


(x/2)2'22  if  x  <  2 
(x-  l)111  ifx  >  2  ’ 


f  2.22  -  V4.94-  1.39e  if  e  <  3 
1  0.81e  —  1.06  if3  <  e  <  18  ’ 


(4)  Degen’s  model: 


V{x) 


(x/2)2'22  if  x  <  2 
(x-  l)111  ifx  >  2  ’ 


J  3.67- V13.44-3.33e  if  e  <  2.65 
1  0.77e  -  0.53  if2.65  <  e  <  18  ’ 


(13) 

(14) 


(15) 

(16) 


(5)  Chang’s  and  Crieri’s  models: 

vh  =  iJ.e2:  (17) 

(6)  Haldar-Miller  model: 

f4.541x  +  0.124  if  0.0388  <  x  <  0.5234 
^(x)  =  2.242x  +  1 .327  if  0.5234  <  x  <  0.746  , 

14.51x  -  7.819  if  0.746  <  x  <  1 .986 


<y(e) 


2.22  —  V4.94  —  1 ,39e  if  e  <  3 
0.81e  —  1.06  if3  <  e  <  18  ’ 


(7)  Hughes’  model: 
ty(x)  =  x,  co^e) 


0.28e  ife<2.52 

0.63e  —  0.89  if  e  >  2.52  ' 


(18) 


(19) 


e 


e 


Fig.  1.  (a  and  b)  Normalized  squared  BLV  as  a  function  of  normalized  thickness  of 
the  shield  calculated  using  different  models:  1:  BRL;  2:  ACE;  3:  NDRC;  4:  Degen;  5: 
Chang  and  Crieri;  6:  Haldar-Miller;  7:  Hughes. 


Fig.  la  and  b  shows  normalized  squared  ballistic  limit  velocity, 
v\{  =  v*  (e)/i^j (5)  =  5)),  as  a  function  of  the  dimen¬ 

sionless  thickness  of  the  shield,  e,  calculated  using  Eqs.  (11 )— (19). 
It  must  be  noted  that  these  models  are  devised  without  requiring 
the  smoothness  of  the  approximating  functions  in  the  endpoints  of 
the  neighboring  intervals.  Consequently,  it  can  be  concluded  that, 
in  the  general  case,  the  dependence  v\{  —  D^(e),  is  described  by  an 
increasing  concave  function  which  approaches  a  linear  function  for 
large  e  (e.g.,  in  the  case  of  the  Haldar-Miller  model). 

Hereafter,  we  use  the  model  given  by  Eqs.  (5)  and  (7)  and  assume 
that  the  function  G(z)  has,  generally,  the  following  form: 


C(z) 


C*(z)  ifO<z<z 
Oi\Z  +  OiQ  if z  >  z* 


(20) 


where  G*(z)  is  an  increasing  concave  function  for  0<Z<Z.,G"(Z»)  = 
0,  coefficients  ao  andoq  satisfy  the  conditions  aiz*  +Qf0  =  G*(z*)  and 
a  i  =  c;(z,). 


2.3.  Model  of  a  layered  shield 


We  study  the  ballistic  properties  of  the  multilayered  shields 
under  the  assumption  that  the  layers  are  perforated  sequentially 
and  independently,  i.e.,  either  the  interaction  between  the  layers 
is  neglected  or  it  is  assumed  that  there  are  relatively  wide  air  gaps 
between  the  layers.  Hence  the  residual  velocity  of  the  jth  plate, 
is  equal  to  the  impact  velocity  of  the  (j+  l)th  plate,  for 

j=l  3-N-l.  In  addition,  the  impact  velocity  of  the  shield  Vimp  is 


equal  to  and  the  residual  velocity  of  the  shield  Vres  is  equal 

to  Vm  =  V(N+1)  (the  notation  V.(N+1)  is  introduced  here  for  conve- 
nience).  Consequently,  the  mathematical  model  of  the  layered  plate 
can  be  written  in  the  following  form: 

f(^p)“F(^p1))  =  G(kr(,))’  J  =  (21) 

where 


imp 


y(N+ 1) 
v  imp 


(22) 


Hereafter,  superscript  in  the  brackets  denotes  the  ordinal  num¬ 
ber  of  the  layer. 


3.  Comparison  of  the  monolithic  and  layered  shields 

In  this  section,  we  compare  the  magnitudes  of  the  BLVs  of  the 
layered  shield  and  of  the  monolithic  shield  having  the  same  thick¬ 
ness. 

Eq.  (21 )  written  for  j  =  1,  2, . . .,  N  and  Eq.  (22)  yield: 

N 

F(Vimp)  -  F(Vres)  =  ^C(fcTffl). 

j=t 


(23) 
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Equation  for  the  BLV  of  the  layered  shield,  Vbb  is  obtained  from 
Eq.  (23)  after  setting  Vimp  =  Vbi  and  Vres  =  0: 

N 

FNbl)  =  (24) 

j=  1 

The  latter  equation  implies  that  in  the  framework  of  the  consid¬ 
ered  model,  the  rearranging  of  the  plates  in  a  multilayered  shield 
does  not  affect  the  BLV. 

In  the  case  of  the  monolithic  shield  with  the  thickness 

Tsum  =  T<1)  +  T<1)  +  ...  +  T<'v)>  (25) 

the  relationship  between  its  BLV,  %,  and  the  thickness  of  the  shield 
can  be  obtained  from  Eq.  (24)  setting  Vbi  =  vbh  N  =  1  and  =  TSUm  • 


(31) 


Eqs.  (30)  and  (31 )  imply  that 

N  -feU) 

A  =  [G'(z)-G'(z  +  Sj_x)]dz.  (32) 

j=2  -1° 


Examining  Eq.  (32)  allows  us  to  arrive  at  the  following  conclu¬ 
sion.  If  a  non  decreasing  in  the  interval  [0,  +  b(2)  h - b 

function  G'(z)  increases  in  some  sub-interval  of  this  interval,  then 
A  <  0,  and  the  monolithic  shield  offers  an  advantage  over  any  lay¬ 
ered  shield.  In  other  words,  this  occurs  if  a  non-convex  function 
G(z)  is  concave  in  some  sub-interval. 


F{vbi)  =  GikTsum).  (26) 

4.  Effect  of  layering  on  ballistic  properties  of  shields 

Taking  into  account  Eq.  (25)  we  can  rewrite  Eqs.  (24)  and  (26)  as 
follows:  4.1.  Formulation  of  the  problem 


N 


F(Vbl)  =  Y)G(bW), 


i= i 


F(vbl)  =  G 


b<f>  =  kT(J\  j  =  1,2, 


(27) 


Consider  the  following  question:  when  layering  does  not  affect 
the  BLV  of  the  shield?  This  occurs  if  F(Vbi)  =  F(%)  for  every  set  b^\ 
. . .,  i.e.,  the  following  equation: 


N 


V'c(f)®)  =  G 

j=l 


(28) 


According  to  the  above  analysis,  ballistic  properties  of  the  shield 
may  be  impaired  by  layering.  Let  us  estimate  the  maximum  reduc¬ 
tion  in  the  ballistic  performance  of  the  shield.  To  do  this  end, 
consider  the  following  problem:  how  one  can  split  a  monolithic 
shield  into  plates  in  order  to  obtain  the  minimum  BLV  of  the  layered 
shield.  We  assume  here  that  the  number  of  the  plates  in  the  shield  is 
given.  The  first  equation  in  Eq.  (27)  implies  that  the  problem  reduces 
to  the  minimization  of  the  following  function: 

N 

A(bm,  b(2), . . . ,  b(N})  =  ^G(b®)-  (33) 

j= i 

Minimization  is  performed  under  the  following  constraints: 


must  be  valid  for  arbitrary  set  of  b^  \  b^2\  . . .,  b^N\  The  solution  of 
the  functional  equation  (28)  is  G(z)  =  cz,  where  c  is  a  constant  (see 
Aczel  and  Dhombres,  1989).  Eq.  (6)  implies  that g(z)  =  G'(z)  =  c  in  this 
case. 

Therefore,  the  BLV  of  the  shield  does  not  change  under  any  lay¬ 
ering  if  and  only  if  the  resistance  in  Eq.  (1 )  does  not  depend  on  the 
instantaneous  DOP  (the  model  in  Section  2.1 )  or  the  BLV  of  a  shield 
is  proportional  to  the  squared  thickness  of  the  shield  (the  model  in 
Section  2.2). 

Now  let  us  compare  the  monolithic  and  layered  shields  hav¬ 
ing  the  same  thickness.  Since  F(z)  is  an  increasing  function,  we  can 
consider  the  difference  of  the  functions  for  two  values  of  the  BLV: 


N 

J2b°)  =  bs'"”’  b(i)>°’  i=h--.,N,  (34) 

j= i 

where  the  total  thickness  of  the  plates  bsum  is  given.  In  the  fol¬ 
lowing  we  consider  an  equivalent  formulation  of  this  problem 
(bsum  =  kTsum)-  minimization  of  the  function 


N-l 


A(bm, b(N~1})  =  J2G(bU))  +  G 

i= i 


min. 


(35) 


N 


4  =  F(Vbl)  -  F(vbl)  =  ^G(fcffl)  _  c 


j=  1 


(29) 


rather  than  the  difference  of  the  BLVs,  Vb[  -  vb[. 

Introduce  the  following  notation:  s,  =  +  h(2)  h - b  b^l\ 

After  some  algebra  we  obtain: 


Minimization  is  performed  under  the  following  constraints: 

N-l 

Y)bU)<b sum,  fa(0>0,  i  =  l,...,N  — 1.  (36) 

j=i 

4.2.  The  case  N  =  2 


G 


‘SN 


Gr(z)dz 


N 

G'(z)dz  +  ^ 
j= 2 


G'(z)dz 


N  „bU) 

=  G(b(1))  +  y^  /  Gfz  +  Sj^dz, 

j=2 


(30) 


First  let  us  consider  the  most  simple  and  the  most  practically 
interesting  case  of  two  layers.  In  this  case  Eqs.  (35)  and  (36)  can  be 
written  as  follows: 

A(bm)  =  G(bm)  +  G{bSUm  -  bm)  min,  0  <  b^  <  bsum.  (37) 
Let  us  calculate  the  derivatives: 

A’(b^)  =  G'(bm)  -  C'(bSum  -  bm), 

A"(b (1>)  =  G"(b(1))  +  G'Xbsum  -  b(1))- 


(38) 
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If  the  interval  with  the  constant  value  of  G\h)  =  G*(z*)  does  not 
exist  (z*  =  oo),  then 


associated  with  the  convexity  of  the  function  A,  and  they  depend 
on  the  signs  of  the  roots  of  the  characteristic  equation: 


A'(bw)  =  c;(b(1))  -  c;(bsllm  -  b(1)), 

A"(b (1))  =  C"(b(1))  +  G"(bSUm  -  b(1>).  (39) 


Since  G*(z)  is  an  increasing  function  for  z>0,  the  con¬ 
dition  for  the  extremum,  A\b^)  =  0,  can  be  satisfied  only 
when  b^  =  bSUm-b^\  i.e.,  equation  A'(b^)  =  0  has  a  single  root, 
=  0.5bSUm-  The  latter  value  provides  minimum  of  A'{b^)  since 
A''(0.5bSUm)  =  2G"(0.5bsum)  >  0.  Therefore,  the  minimum  BLV  is 
attained  when  the  monolithic  plate  is  separated  into  two  plates 
having  equal  thicknesses. 

Assume  now  that  z*  <  oo,  i.e.,  the  interval  with  the  constant 
value  of  G\z)  =  G*(z*)  does  exist.  Let  us  show  that  the  equation 

can  be  satisfied  only  in  two  cases:  b(1)<z*  and 
fi(2)<Z*  or  b >z*  and  b^  >z*.  The  proof  is  given  by  the  following 
two  sequences  of  inequalities:  the  first  one,  <  z*  =>>  G'(fi(1))  < 
G'(z *)  =>-  G'{b (2))  <  Gr(z*)  =>  fi(2)  <  z*  and  the  second  one,  > 
z*  =>  G'(ib(1))  =  G'{z *)  =>  G'(h(2))  =  G'(z*)  =>  b(2)  >  z*. 

If  bsum  <  2 z*,  then  only  the  first  possibility  <z*  and  fi(2)  <z*) 
can  be  realized  in  the  point  of  the  extremum  of  A(b^).  Conse¬ 
quently,  Eq.  (39)  has  the  form  of  Eq.  (38)  and  the  single  minimum 
of  zl(fi(1))is  attained  for  =  0.5bsum)- 

If  bsum  >  2z*,  then  the  condition  of  extremum,  A'(b^ )  =  0,  implies 
the  inequalities  fi(1)>z*  and  b(2)  =  bsum-b^  >z*.  Since  A"  =  0  for 
such  if/1),  it  is  convenient  to  consider  the  following  relationship  for 
G'{b™): 

G\b M)  =  A\b M) 


G^b™)  -  G'*(z*)  <  0 

0 

G'Sz*)  -  G'Sbsum-b^)  >  0 


if  0  <  <  z* 

ifz*  <  fi(1)  <  bsum-Z* 

if  bSum~Z*  <  b^  <  bSum 


(40) 


Function  A'(b^)  decreases  when  0<b^<z*,  remains  constant 
when  z*  <  -  z*  and  increases  when  bSUm  -  z*  <  )  <  bSUm  )• 

Since  z*  <  0.5  bSUm  <  bsum  -2*,  maximum  of  Ar(fi(1))  equals  to  zl(0.5 

fisum  )• 

Consequently,  we  established  and  proved  the  following  asser¬ 
tion.  If  bsum  <2 z*,  then  the  minimum  BLV  is  attained  only  when 
the  monolithic  plate  is  separated  into  two  plates  having  equal 
thicknesses.  If  bSUm  <  2z*  then  the  minimum  BLV  is  attained  for 
all  z*  <  <  bsum  and  z*<b(2)  <bSUm  including  the  particular  case 

of  splitting  the  shield  into  two  plates  with  equal  thicknesses. 


4.3.  General  case 


Let  us  investigate  the  case  with  arbitrary  number  of  layers.  Since 
the  first  and  the  second  order  derivatives  of  G  exist,  zl  is  a  suffi¬ 
ciently  smooth  function.  The  equations  for  the  first  and  the  second 
order  derivatives  of  A  read: 


H{v)  =  99bfv)  =  G'(b^)  G'(bW), 

(41) 

r)2  A 

dbMdbM  =i"C"(fcM)  +  G"(bM)’ 

(42) 

where  v  =  l,...,N-l,/z  =  l,...,N-l, 

Sv/X  =  {  0  if  V  +  M  ’  biN)  =  bsum  ~  YlbU)- 

k  j= i 

(43) 

The  Hessian  (the  matrix  with  the  elements  plays  the  role 

of  the  second  derivative  of  A,  considered  as  a  function  of  many  vari¬ 
ables  (see,  e.g.,  Korn  and  Korn,  1968).  Properties  of  the  Hessian  are 


H(T2)  ...  hV>n~  1} 

H(2,1)  h(2,2)_^  ...  Ht2^-1) 

•  ••  •••  •••••• 

H(N-1,1)  H(N- 1,2)  _ 


(44) 


Since  the  last  term  in  Eq.  (42)  is  the  same  for  all  H^v^\  this  term 
can  be  dropped  out  in  the  determinant  in  Eq.  (44).  Nonzero  ele¬ 
ments  remain  only  at  the  diagonal,  and  the  characteristic  equation 
reads: 


N-l 

r[[G"(fa(i'))-??]=o. 

i= 1 


(45) 


The  domain  determined  by  the  linear  inequalities  (36)  is  convex. 
Since  G"(z)  >  0  in  this  domain,  all  the  roots  of  Eq.  (45)  are  nonneg¬ 
ative,  and  hence  the  Hessian  is  positive-semidefinite.  In  classical 
mathematical  analysis  such  functions  are  known  as  non-convex 
functions. 

Let  us  prove  that  the  considered  problem  always  has  a  unique 
solution  if 

bsum  <  Nz*.  (46) 

Using  Eq.  (41),  let  us  rewrite  the  conditions  =  0(v  = 
1 , . . . ,  N  -  1 )  in  the  following  form: 

C'(b(1>)  =  C'(b(2>)  =  •  ■  •  =  G'(b(N~V)  =  G'(bm).  (47) 

Eq.  (46)  implies  that  there  exists  at  least  one  b ^  <  z*.  Taking 
into  account  the  property  of  the  function  G(z)  we  conclude,  similar 
to  the  case  N  =  2,  that  all  <  z*  and  Eq.  (47)  have  unique  solution: 

bd)  =  6<2>  = . . .  =  gw- 1)  =  m  =  .  (48) 

N 

Since  all  <  z*,  then  G"(fi(v))  >  0,  Eq.  (45)  has  only  positive 
roots  and  the  Hessian  is  positive-definite.  Therefore,  A  is  strictly 
concave  at  the  point  b and  this  point  is  the  point  of  a 
local  minimum.  In  convex  programming  (the  considered  problem 
belongs  to  the  class  of  problems  considered  in  convex  program¬ 
ming),  a  local  minimum  (if  it  does  exist)  is  also  a  global  minimum 
(see,  e.g.,  Rockafellar,  1970).  Therefore,  it  can  be  concluded  that  the 
global  minimum  of  function  A  is  attained  at  the  point  determined 
by  Eq.  (48),  i.e.,  when  the  thicknesses  of  the  plates  are  equal. 

In  the  case  when 


bsum  >  Nz*, 


(49) 


there  exists  a  non-empty  sub-domain  of  the  domain  determined  by 
Eq.  (36)  that  is  associated  with  values  b for  which  the  BLV  is  min¬ 
imum.  This  sub-domain  is  determined  by  the  following  inequities: 


N-l 


Z*  < 


b(v)  <  b 


—  “sum, 


^  ^b^  <  bsum  —  Z* ,  V  —  1 ,  .  .  .  ,  N  —  1 .  (50) 


j= 1 


The  structure  of  this  sub-domain  for  N=3  is  showed  in  Fig.  2. 

Fig.  3  shows  maximal  decrease  of  the  BLV  as  a  function  of  the 
number  of  layers  for  the  NDRC  model.  In  order  to  quantify  the 
decrease  of  the  BLV  we  use  the  normalized  BLV,  %(N) /%( 1 ),  where 
%(N)  is  the  BLV  of  the  layered  shield  with  N  equal  layers  and  %( 1 ) 
is  the  BLV  of  the  monolithic  shield  with  the  same  thickness.  When 
v2hl  is  proportional  to  en  ( n  =  1.5  for  the  BRL  model,  n  =  2  the  Chang’s 
model  and  Crieri’s  model),  the  following  simple  relationship  can  be 
obtained: 


VblW 
vblC 1) 


=  N(l-n)/2i 


(51) 
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Fig.  2.  Sub-domain  of  the  thicknesses  of  the  plates  for  a  three-plate  shield  with 
values  of  and  that  provide  the  largest  reduction  of  BLV  by  layering. 


N 

Fig.  3.  Decrease  of  the  BLV  with  increasing  the  number  of  layers  for  the  NDRC  model. 

Eq.  (51)  and  Fig.  3  show  that  layering  can  cause  a  considerable 
reduction  of  the  ballistic  performance  of  a  concrete  shield.  Clearly, 
the  quantitative  results  obtained  using  the  approximate  models  can 
be  of  interest  primarily  as  a  basis  for  qualitative  conclusions. 

5.  Concluding  remarks 

Using  a  wide  class  of  semi-empirical  models,  we  investigated 
the  effect  of  layering  of  concrete  shields  and  established  a  number 
of  general  assertions  concerning  ballistic  properties  of  the  shields. 
Omitting  the  technical  details,  the  main  conclusions  of  this  study 
can  be  formulated  as  follows: 

( 1 )  The  BLV  of  the  multilayered  shield  does  not  depend  on  the  order 
of  the  plates  in  the  shield. 

(2)  Monolithic  shield  performs  better  than  any  layered  shield  with 
the  same  thickness. 

(3 )  The  largest  decrease  of  the  BLV  occurs  when  the  shield  is  divided 
into  the  plates  having  the  same  thickness. 


Note  that  we  did  not  specify  a  particular  relationship  between 
the  BLV  and  the  thickness  of  the  shield,  but  used  only  some  general 
properties  of  this  dependence. 

Clearly,  the  obtained  theoretical  predictions  require  experimen¬ 
tal  validation. 
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